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Abstract 

In this paper we consider in detail the connection between the problem of a 
polymer in a random medium and that of a quantum particle in a random 
potential. We are interested in a system of finite volume where the polymer 
is known to be localized inside a low minimum of the potential. We show 
how the end-to-end distance of a polymer which is free to move can be ob- 
tained from the density of states of the quantum particle using extreme value 
statistics. We give a physical interpretation to the recently discovered one- 
step replica-symmetry-breaking solution for the polymer (Phys. Rev. E61, 
1729 (2000)) in terms of the statistics of localized tail states. Numerical solu- 
tions of the variational equations for chains of different length are performed 
and compared with quenched averages computed directly by using the eigen- 
functions and eigenenergies of the Schrodinger equation for a particle in a 
one-dimensional random potential. The quantities investigated are the radius 
of gyration of a free gaussian chain, its mean square distance from the origin 
and the end-to-end distance of a tethered chain. The probability distribution 
for the position of the chain is also investigated. The glassiness of the system 
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is explained and is estimated from the variance of the measured quantities. 
PACS number (s): 36.20.Ey, 05.40-a, 75.10.Nr, 64.60.Cn 
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I. INTRODUCTION 



The behavior of polymer chains in random media is a well studied problem [1-8] that 
has applications in diverse fields. Besides the polymers themselves this problem is directly 
related to the statistical mechanics of a quantum particle in a random potential |J, the 
behavior of flux lines in superconductors in the presence of columnar defects [II], IT], and the 



problem of diffusion in a random catalytic environment H. It was found in Ref. [0-^|,|7J that a 
very long Gaussian chain, immersed in a random medium with very short range correlations 
of the disorder, will typically curl up in some small region of low potential energy. The 
polymer chain is said to be localized and for long chains the end-to-end distance becomes 
independent of chain length (R 2 ~ L°). For short chains the end-to-end distance scales 
diffusively (R 2 ~ L). In a related paper ||, it was found, using the replica approach, that a 
quantum particle exhibits glassy behavior at low temperatures. The low temperature limit 
for a quantum particle translates into the long chain limit for polymers. It implies that 
the free energy landscape of a long chain is typically very complicated and possesses many 
metastable states. 

Recently a new variational solution has been found for gaussian chain embedded in a 
random 5-correlated potential |7]]. This solution, which involve replica-symmetry-breaking 
(RSB) at the 1-step level, gives rise to the correct behavior of the end-to-end distance of the 
chain as predicted by the heuristic free-energy estimates of Cates and Ball [[| . It predicts the 
subtle log V dependence of the chain size, which was not present in the variational solution 
of Edwards and Muthukumar |2| due to the fact that their solution did not contain enough 
variational parameters and hence did not reflect the translational invariance of the original 
hamiltonian. 

Our aim in this paper is two fold. The first goal is to strengthen and elucidate the 
connection between the localization of a polymer in a random medium and that of a quantum 
particle. There is an extensive literature on electron localization that we would like to build 
bridges into. The second goal is to give further interpretation and a better physical picture 
of the recently discovered variational solution for a polymer in a random medium W\. It 



turns out that these two goals intertwine together by using the properties of the solutions 
of the Schrodinger equation for a quantum particle in a random environment. 

Using the path integral mapping between the partition function of a Gaussian polymer 
chain and the imaginary time Schrodinger equation, we show that the glassy behavior of 
the polymer chain can be understood by studying the eigenf unctions of the Schrodinger 
equation with a random potential. In particular, the phenomenon of Anderson Localization 
ll"2l , |l3H is crucial to understanding the glassy phase. We present evidence that shows that 
polymer localization and glassy behavior can be traced to the dominance of exponentially 
localized tail states. We also explore various connections between typical glassy behavior, 
such as non-ergodicity and the existence of many metastable states, to the properties of the 
eigenfunctions of the Schrodinger equation. 

In order to describe the glass phase analytically we utilize the variational replica approach 
described in f7j and ||. We further investigate and interpret the transition discovered there 
between a replica symmetric phase to a 1-st ep replica symmetry breaking phase at a critical 
chain length. The 1-step breaking solution describes the localized (glassy) phase of the 
polymer and corresponds to the dominance of localized tail states. In the long chain limit 
the lowest tail state in each realization of the disorder will dominate the partition sum. 
On the other hand, the replica symmetric phase corresponds to the case when a multitude 
of extended states dominate the partition function. Our results are substantiated by a 
numerical solution of the problem by solving the imaginary time Schrodinger equation on a 
lattice in one spatial dimension. 



II. THE MODEL 



The simplest model of a polymer chain in random media is a Gaussian (flexible) chain 



14j| in a medium of fixed random obstacles ||. In this paper we do not include a self-avoiding 



interaction. This model can be described by the Hamiltonian 



H 



du 



M fdR{u) 
Y l du 



+ |R» + V(R(u)) 



(2.1) 
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where R(it) is the d dimensional position vector of a point on the polymer at arc-length 
u (0 < u < L), and where L is the contour length of the chain (number of links). The 
random medium is described by a random potential l^(R) that is taken from a Gaussian 
distribution that satisfies 

(V(R)> = 0, (V(R)V(R'))=f((R-R>) 2 ). (2-2) 

The harmonic term in the Hamiltonian is included to mimic the effects of finite volume. 
This is important to ensure that the model is well defined, since it turns out that certain 
equilibrium properties of the polymer diverge in the infinite volume limit (fi — > 0). The 
function / characterizes the correlations of the random potential, and will depend on the 
particular problem at hand. The parameter M is equal to d/(/3b 2 ), where f3 = (fcgT) -1 , and 
where b is the Khun bond step. 

In this paper we will consider Gaussian distributed random media defined by the corre- 
lation function 

/ (( R " R ') 2 )) = exp (— (R - R') 2 /e 2 ) . (2.3) 

Here g determines the strength of the disorder and the parameter £ controls the correlation 
range of the random media. For £ large (long range correlations) we expect to recover the 
exact results presented in Ref . || . When £ is small then we have the more physically relevant 
case of short range correlated media. In the limit of £ — > 0, / approaches a <i-dimensional 
5-function. This was the case studied in Refs. . In the present paper we investigate only 
the case of short-range correlations, i.e. the case of small £. 

Once we have defined the Hamiltonian for any chain configuration R(m), we can write 
the partition sum (Green's function) for the set of paths of length L that go from R to R' 
as 

/•R(L)=R' 

Z(R,R;L)=/ [dR(u)]exp(-pH). (2.4) 

</R(0)=R 

All the statistical properties of the polymer will depend on the partition sum. For 
instance, the end-to-end distance (or radius of gyration) of a polymer chain that is free to 
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move is given by 



/ dRdR'(R - R) 2 Z(R, R; L) 



< R ^» = { T^izikM) ' ) • (2 ' 5) 

where the overbar stands for the average of the ratio over the realizations of the random 
potential. This average is referred to as a quenched average, as opposed to an annealed av- 
erage, where the numerator and denominator are averaged independently. In some previous 
studies it has been customary to replace the quenched average of certain quantities (and in 
particular Rp) by the more analytically tractable annealed average. However, this replace- 
ment can be justified in the case of Eq. (|2.5|) only when the system size is strictly infinite, 
since only in that limit can the polymer sample all of space and find the most favorable 
potential which will be similar to its environment in the annealed case. The main problem 
with this approach is that in practice we always deal with finite size systems, and in general 
it is difficult to assess how big the system size has to be so that the annealed average is a 
good approximation to the quenched average. In addition, the time it takes the chain to 
sample a large volume is exceedingly long and unreachable over a reasonable experimental 
time. In this paper we will always deal with an explicitly finite system (both analytically 
and numerically) and so we will only compute quenched averages. For very large but finite 
systems the free chain is said to be localized in the sense that its configuration space is 
dominated by a single configuration in which the chain is being trapped in a single small 
neighborhood. The size of the chain in an uncorrelated random potential is given by 0,0] 

R F oc (g\nV)- 1/(d ^\ (2.6) 

where V is the volume of the system. The depth of the well entrapping the chain is approx- 
imately 

K™„~G?lnV) 2/(d - 4) . (2.7) 
This is the binding energy per monomer. Thus the binding energy of the chain is given by 
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Another quantity of interest is the averaged mean squared displacement of the far end 
of a polymer with one end that is fixed at the origin. This is a measure of the wandering of 
a tethered polymer immersed in a random medium. This quantity can be written as 



/ dRR 2 Z(0,R; L) 



(RfW) = [ J dRZ(0, R;L) )' (2 ' 8) 
It is important to distinguish between these two quantities since they measure different 
properties of the polymer chain. Cates and Ball || and Nattermann |3J discussed the 
different behaviors of these quantities for a short-range correlated random potential using 
heuristic arguments. In a recent paper || we have found analytically that the scaling 
properties of Rp and Rt are also very different when the disorder has long range quadratic 
correlations. This is because in the case of a tethered chain, the free end of the chain seeks 
out favorable regions in the random medium that are typically very far from the origin. 
When both ends are free then the entire chain will simply curl up in a favorable region in 
the random medium, and the end-to-end distance will not scale as fast. For the tethered 
chain the quenched and annealed averages are not expected to coincide even in the infinite 
volume limit, since a tethered chain can only sample a finite volume of size L d , and even for 
large L the environment near its tail can not be adjusted since it is immobile. 
Yet another quantity of interest is 



/ dKR 2 Z(K, R; L) 



(RqW) = \ J dRZ(R, R;L) )' ^ 
which measures the distance from the origin (of the harmonic potential) to the average 
(center of mass) position of a free chain (more precisely of a free loop). It is also relevant to 
the related problem of a quantum particle in a random potential P,|T5| . The relation between 
the polymer and the quantum particle problems arises from the fact that the partition sum 
of a polymer chain can be mapped to the density matrix of a quantum particle. The mapping 
HH is given by 

P^l/h, L^(3h. (2.10) 

Then p(R, R'; (3) = Z(R, R'; L = /3h, (3 = 1/ti) is the density matrix of a quantum particle at 
inverse temperature j3. Note that the variable u is now interpreted as the Trotter (imaginary) 



time, and M as the mass of the quantum particle. Under this mapping (Rq(L)) can be 
interpreted as the average mean squared displacement from the origin of a quantum particle 
in a random plus harmonic potential centered at the origin. 

Finally we would like to remind the reader the definition of the density of states for 
a quantum particle in a random potential. We first define the quantity p((3) which is the 
two-sided Laplace transform of the density of states p(E): 

/oo 
exp(-/3E) p(E) dE. (2.11) 
-oo 

This function is given by 

ptf) = lim i I dRZ(R,R;£). (2.12) 

V^oo y J 

We will make use of these quantities in the following sections. 

III. THE PATH INTEGRAL MAPPING 



The partition sum of the polymer chain (|2.4j) can be mapped to an imaginary time 
Schrodinger equation. This mapping (see Ref. |1J Eqs. (3.12)-(3.18)) is given by 

»R(L)=R 
'R(0)=R' 

where 



l>R(L)=R 

Z(R,R';L) = [dR(u)]exp(-pH{R(u)]) = (R\exp(-/3LH)\R'), (3.1) 

Jn(o)=R' 



H= l —^- + ^-R 2 + V{R). (3.2) 

So for a given realization of the random potential the polymer partition sum can be expressed 
as a matrix element of the imaginary-time evolution operator. The matrix elements can be 
expanded in a complete set of eigenfunctions of the Hamiltonian operator to yield 

oo 

(R|exp(-/3L#)|R> = ^exp(-/5LE m )$ m (R)$;(R), (3.3) 

m=0 

where 

£$ m (R) = E m $ m (R). (3.4) 
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It is clear from the above relation that all properties of the polymer chain will depend on 
the eigenvalues and the eigenfunctions of the Hamiltonian operator H. For the problem at 
hand we will have to understand these properties for the case when the potential V(R) is 
random and with the correlations given in Eq. ( |2.3| ) with a small £. 

The Schrodinger equation with a random potential is a well known problem that has 



been intensely studied for a long time |12| , [13| , |17| , P5| . The main property that we will use 
is that when V(R) has short range correlations (i.e. correlation length is shorter than any 
other length scale in the problem), and if the system size is infinite, then in any dimension 
all eigenstates with energy bellow a critical energy Em ( refered to as the mobility edge) are 
exponentially localized in the form 

$ m (R) -expHR-Rj/40. (3.5) 

Here R m is the localization center of the m th state, and i m is the localization length of that 
state. The localization length satisfies 

l/£ m = (3^2M\E m \, (3.6) 

for E m «C 0, i.e. deep in the tail region. Intuitively, it is easy to verify these last two relations 
for the solution of the Schrodinger equation for a particle in a one dimensional, non-random, 
attractive 5-function potential. This can be thought to represent a local minimum of the 
random potential. In <i-dimensions one can similarly consider a potential of the form — vS(r) 
and the lowest energy solution of the radial equation (with zero angular momentum) also 
satisfies these relations. For E > Em extended states exist when d > 2. For d = 1, 2 there is 
no mobility edge and all states are exponentially localized. The states with energies E > Em 
are called extended since they are no longer localized but are spread over a finite fraction of 
the system. Also, it is known that the eigenvalues of the localized states are discrete, while 
the eigenvalues of the extended states form a continuum. 

For finite system size, or if /i ^ in the Hamiltonian given in Eq. ( p.2|) , the above 
discussion has to be modified. First, the eigenfunctions are always discrete in any dimensions. 
But even in one dimension as the energy increases the width of the localized states eventually 
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becomes comparable to the system size and thus a localized particle of that energy can 
go from one end of the sample to the other. Thus the distinction between localized and 
extended states becomes blurred for a finite system at energies much above the ground 
state. Nevertheless, there will still be a qualitative difference between the low energy tail 
states and the higher energy states with large localization lengths. For a finite system (or 
when /i ^ 0) and for any given realization of the random potential there will always be a 
lowest energy state which is by definition the ground state for that realization. If the volume 
is large (or \x small) it will correspond to one of the deep tail states of the infinite system. 

In order to study the effect of the eigenfunctions on the physical properties of the polymer 
chain, we simply apply Eq. Q3.3|) . It is clear that when L (or f3) is large the partition sum 
is dominated by contributions due to the low energy localized states, while if L is small 
then most of the contributions will come from the continuum of extended states. As L is 
slowly increased then fewer and fewer localized states contribute to the sum, until we attain 
ground state dominance at some very large L. For simplicity let us concentrate on the 
diagonal elements of the evolution operator. Then 

Z(R,R;L) = ^exp(-/3L£ m )|$ m (R)| 2 , (3.7) 

m 

which is proportional to the probability of finding a loop of length L passing through the 
point R. Since the eigenstates for low energies are localized then for large L we can write 

Z(R,R;L) ~^exp(-/3L£ m )exp(-2|R-R m |/£ m ). (3.8) 

m 

This implies that the probability of finding a long loop passing through R is concentrated 
around the localization centers of the low energy eigenstates. Consequently, a very long 
polymer chain will most likely be found at the ground state localization center. When L is 
small then the partition sum will be dominated by the extended states, and so the probability 
of finding a loop at R should be fairly uniform throughout the system. In a later section we 
will analyze in detail the evolution of the partition sum with L by solving the Schrodinger 
equation on a lattice in d = 1. 
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All the physical properties of the polymer chain can be expressed in terms of the eigen- 
states of Schrodinger equation. For instance we can write the end-to-end distance for a given 
realization of the random potential as 

2 E m L/rfRR 2 ^(R) - |J rfRR$ m (R)| 2 ) exp(-(3LE m ) 

where a m = J dR$ m (R), and where (-)y refers to a configurational average for the case of 
a fixed realization of random potential. When L is large enough so that (E\ — E gs )L >> 1, 
where Ei is the eigenvalue of the first excited state, then only the ground state contributes. 
In this case we have 

, 2 . r .. _JdRR^ gs (R) //dBM^R) 



< r ^- 2 t^w-Ht^^J • (310) 

where <& ffS (R) is the ground state eigenfunction. It can be shown that the ground state wave 
function is positive definite and so in the large L limit (R^(L)) 1 ' 2 can be interpreted as the 
width of the ground state eigenfunction. Assuming the ground state has the form given in 
Eq. ( p.5|) , we can write (R F (L)) V = 2d(d+ l)£g S , where £ gs = £q is the localization length of 
the ground state. Upon averaging over all realizations of the random potential we get that 



(R F (L)) = 2d(d + l)£g S , and so the quenched average of the end-to-end distance, in the long 
chain limit, is proportional to the square of the average localization length of the ground 
state eigenfunction. We see that for \i — > the average is taken over the tail states of the 
Schrodinger equation. 

When one end of the polymer is tethered to the origin then the end-to-end distance can 
be expanded in eigenfunctions to yield 

, B3m v _ E m IdRR 2 ^ m (R)^ m (o)eM-PLE m ) 

{tlAL))v ~ E m /^(R)$ m (0)exp(-/3LE m ) • [6AL) 
The presence of the term $ m (0) is crucial since if the eigenstates are localized then $ m (0) ~ 
exp(— |R m |/£ m ), which means that eigenstates localized far away from the origin may not 
contribute to the sum even if they may have very low energies. When L is not very large 
the extended states will dominate the sums in Eq. ( p.ll|) and the behavior of the end- 



to-end distance should be diffusive ((Rf.(L))y ~ L). When L is large, but before the 
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onset of ground state dominance, the sums in Eq. ( ft.lip will be dominated by localized 



eigenstates that are centered close to the origin. The free end of the polymer chain will 
hop between localization centers as L increases. This behavior has been investigated using 



Flory arguments in Ref. and both authors found a weakly sub-ballistic ( (R|,(L)) ~ 
L 2 /log(L) 7 ) behavior of the mean squared displacement. For a finite system and when L is 
large enough the ground state will always dominate and we get (R^(L))y = d{d+l)£ 2 gs +~R? gs . 
Since the distance from the ground state to the origin will typically be much larger than the 



localization length, we find upon averaging that (R^(L)) rs R 2 S , and so the quenched average 
of the end-to-end distance simply converges to the distance to the localization center of the 
ground state. This distance is of the order of the size of the system. Thus the interesting 



sub-ballistic L-dependence of (R^(L)) arises from the contribution of excited states and not 



from ground state dominance as for the case of (R? F (L)). 

The mean squared displacement defined in Eq. (|2.9|) can also be expanded in eigenfunc- 
tions as 

/W 2 mN E m / rfRR 2 |<MR)| 2 exp(-(3LE m ) 
In the large L limit we find that (Hq(L)) v = d(d + l)£g S /4 + R 2 s . This is very similar to the 



case discussed in the previous paragraph and we find upon averaging that (Rq(L)) « R : 



gs ■ 



Which implies that (Rq(L)) becomes independent of L. When L is small then the extended 
states dominate and we expect a behavior as if there is no random medium. For a polymer 



confined by a quadratic potential we find that (Rq(L)) = d(PfiL)^ 1 . 

The relationship between the eigenstates of the Schrodinger equation and the partition 
sum of the polymer chain can also be used to understand the behavior of the sample-to- 
sample variations of various physical properties of the chain. The sample-to-sample varia- 
tions are important in order to assess whether the average (over many samples) of a physical 
quantity is an accurate measure of that quantity for a typical sample. For instance, the 



sample to sample variation of (Hp(L)) v is defined as A F = (R F (L)) — (R|,(L)) . In the 



same way we can define the sample-to-sample variation = (R^(L)) 2 — (R^(L)) and 



Aq = (Rq(L)) — (Rq(L)) . So if the sample-to-sample variation of a physical quantity is 
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large then that implies that the probability distribution for that quantity is highly skewed 
and strongly dependent on the random sample. Systems which exhibit such behavior are 
said to have glassy properties. For simplicity if we assume ground state dominance (L very 
large) then A^ is just a measure of the sample-to-sample variation of the localization length 
of the ground state. Also and Aq are both a measure of the sample-to-sample variation 
of the distance from the origin to the ground state localization center. The crucial point 
here is that the sample-to-sample variations are dependent on which eigenstates dominate 
the partition sum, and so there should be a connection between glassy properties and the 
behavior of the eigenstates. 



IV. THE REPLICA VARIATIONAL APPROACH 

In this section we review the replica approach that was used in Refs. ]7| and || to 
compute quenched averages of various physical properties of the polymer chain. Our goal 
is to give an interpretation for this formalism in terms of the localized states picture of the 
corresponding Schrodinger equation. 

In order to compute the quenched average over the random potential we apply the replica 
method. We first introduce n-copies of the system and average over the random potential 
to get 



where 



Z n ({Ra}, {R'J; L) = Z(R U R£; L) ■ ■ ■ Z(K n , R' n , L) 

r.R a (L)=R4 n 

l[[dR a }exp(-(3H n ), 

R a (()): R a a=l 



-£ I du [ L d U 'J2f((Ra(u)-R b (u')) 2 ). 

1 Jo Jo nh 



(4.1) 



(4.2) 



The averaged equilibrium properties of the polymer can now be written in terms of the 
replicated partition sum Z n ({R a }, {R a }; L). For instance, the mean squared end-to-end 
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distance denned in Eq. Q2.5Q can be written in as 



/ n dR a n d K (Ri - Ri) z n ({R a }, {R'J; L) 



(R|(L)) "-0 /n^an^'^n({Ra},{R' a };L) a ' ' ^ 

Thus, the averaged equilibrium properties of the polymer can be extracted from an n-body 
problem by taking the n — > limit at the end. This limit has to be taken with care, by solving 
the problem analytically for general n, before taking the limit of n — > 0. Unfortunately the 
replicated partition sum cannot be evaluated analytically and a variational approach has 
been used in Refs. |7]J| to make further progress. The procedure is to follow the work of 



Feynman |L6| and others |T9| , |20[| and model H n by a solvable trial Hamiltonian h n which is 



determined by the stationarity of the variational free energy 

n(F) = (H n - h n ) hn -\^J MRi] • ■ ■ [dRn] exp(-^ n ). (4.4) 

Note that the variational free energy also depends on the boundary conditions on the 
polymer chain. If we are interested in the case when one end is fixed then the parti- 
tion sum should be over paths with one end fixed. In this case the path integrals should 
be evaluated using / dRiJ^J^'idRi]. When both ends are free one should use in- 
stead J dRjfiR- f^lfy Jft 4 [dRi] ■ Alternatively, one can evaluate the partition sum with 
J dRi J^tyJ^ 1 [dRi] , which yields the free energy of a polymer loop of contour length L. 
The parameters that minimize the variational free energy for this case should be the same as 
for the case when both ends are free since in both cases the polymer loop can be anywhere 
in the system. These are the boundary conditions that were used in Ref. 0, where only the 
polymer chain that is free to move was investigated. 

The quadratic trial Hamiltonian has been parametrized J/J by 



~77 I du I du'^pab (R a (u) - Rb{u')) 
4jL Jo 



(■1.0) 

ab 

where the matrix elements p a b are the variational parameters. The physical motivation for 
this ansatz is that the replica-replica interaction in the original Hamiltonian is modeled by 
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a quadratic interaction which can be different for different replica pairs. Also, as noticed in 
0, the quadratic interaction has the same translational invariance as the interaction term 
in the original Hamiltonian. The case when the replica-replica interactions are the same 
(Pab = constant) corresponds to the case of long range quadratic correlations which was 
solved exactly in ||. If we expand and simplify the quadratic interaction we can rewrite the 
trial Hamiltonian as 

+^7 I du I ^TpAM-RjM, (4.6) 

lL J ° J ° ab 

where A = [A — J2bPib- Here, A is assumed to be independent of the replica index a, as is 
the case if p is a hierarchical matrix. In Ref. A was treated as an independent variational 
parameter, and the condition A = \i — 'Yl lb p a b emerged automatically as a result of the 
translational invariance of the disorder dependent term in the original hamiltonian H n . 
In the 1-step RSB scheme, the matrix p ab can be parametrized as (p,p(x)) with 

{—s < x < x c 
(4.7) 
— Si X c < X < 1 , 

and where x is Parisi's replica index. In J7] p has been denoted by A — Ai. Thus, five 
variational parameters were used: A, Ai, sq, s± and x c . The variational free energy was 
expressed as a function of these variational parameters. Taking the partial derivative with 
respect to the variational parameters and equating them to zero one gets five non-linear 
equations which could be solved analytically when L was large and \i was small 0. In 
this paper we started from the same free energy and took its partial derivatives without 
simplifying the expressions for large L. As mentioned above because of the consequences of 
translational invariance we could reduce the number of parameters and equations to four, 
since 

Ai + x c s + (1 - x c )si = fi . (4.8) 
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We proceeded to solve the equations of stationarity numerically using a standard iterative 
method P^fl . Once the variational parameters have been found we can obtain expressions 
for various physical quantities. For instance, we find that 

2d 



{RUL)) = TTml — 1 7a \ v (4 ' 9) 

and also that 

d ( 1 1 s (1 — x c )S \ d (/XL 



< R ^» = + ii&TTj) + wm coth \ \l ~m~2 J ' < 4 - 10 ' 

where we have put £ = x c (s\ — s ). Details of the numerical results will be given in the 
next section. For large L and small \x it can be seen using the results obtained in ]7| that to 
leading order 

2d 



(Rl(L)) = I_ , (4.11) 



and 



with 



and 



^4/(4 d) (4+d)/(4-d) / |, i\4/(4-d) 1Q x 

A = ( 27r )2d/(4-d) ^ M ) IM) > ( 4 - 13 ) 



i / Jd-2 \ -l/(4-d) 



Xc = Z vfivrp l 1 * 1 /^ ) ■ (LID 



V. NUMERICAL PROCEDURE 

We check the validity of the analytic solution by numerically computing the quenched 
average of various physical properties of the polymer. This is computationally intensive be- 
cause all quantities will have to be averaged over many realizations of the random potential. 
In this paper we will only concentrate on the case d — 1. Although this does not correspond 
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to a physical polymer (d = 3) we will still be able to check the validity of our analytical 
results for the special case d — 1. 

We evaluate numerically the right hand side of Eq. ( |3 . 3| ) by solving the Schrodinger 
equation on a one dimensional lattice of N sites ||22|| . The lattice Hamiltonian is an iV x iV 
matrix with matrix elements given by 

= ~2M^W (< ^ +1 + Si+1 ' j) + (2 A2( * " N/2)2 + V{{) ) kj (5 ' 1} 

where the lattice spacing is A = S/N, and where S is the system size. Since we are 
interested in the continuum limit A will be kept small. Note that the index i corresponds 
to the position Ri = Ai. We impose hard wall boundary conditions at the end of the 
lattice. The eigenvalues and eigenvectors can now be found directly by diagonalizing the 
matrix using a standard numerical routine p2| . Once these are known we can construct 
the partition sum at any value of L using Eq. ( |3.3| ). The partition sum can then be used 
to compute the quantity of interest, such as (R F (L)) V . We repeat this procedure a large 



number of times and average the results to get a numerical approximation to (R, F (L)). 

The correlated Gaussian random potential described by Eq. (|2.3|) is modeled by a se- 
quence of N numbers {V^(z)}j = i r . ) jy that obey (V^(i)V^(i + I)) oc exp (— A 2 / 2 /£ 2 ). These 
numbers are then placed on a lattice of iV sites in the given order. To generate such num- 
bers we use an established method for generating correlated random numbers. The details 
of this method are described in Ref. [01. 



VI. RESULTS AND DISCUSSION 

A. Numerical and Analytical results for (R? F (L)) and (Rq(L)) 

Using the method described in the previous section we compute (R F (L)) on a lattice of 
size N = 200 (S=40). Here, we concentrate on the case where the random potential has 
very short range correlations. We generate random potential samples with correlation length 
£ = 1/ y/2. In Fig. [1] we show a typical sample of 200 numbers from such a distribution. 
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FIG. 1. A typical random sample with £ = l/y/2 

Notice that the short range correlation leads to a rugged random potential landscape. The 
average mean squared displacement (Hp(L)) is computed by averaging over 10000 samples, 
and the error is estimated by computing the standard deviation of 10 sets of 1000 samples. 

We were able to find numerical solutions to the non-linear stationarity equations. We 
found that for a given set of parameters there is a chain length L c (which depends on the 
strength of the disorder) such that for < L < L c there is only a replica symmetric solution. 
This is the case when the variational parameters satisfy x c — 1 and Sq = Sx- For L > L c 
there is still a replica symmetric solution but we also find an additional replica symmetry 
breaking solution. So in this regime we find an additional solution such that < x c < 1 and 
sq 7^ sx- In order to decide which solution correctly describes the physics in that regime we 
compare their respective predictions to the lattice computation of (R F (L)) and (Rq(L)). 

In Fig. |2| we plot the mean squared displacement (RJ,(L)) vs. L for a given set of 
parameters. We plot this quantity using the lattice result, and also using the two predictions 
of the variational method. Note that in the labels of the plots the average over the disorder 
is denoted by a second set of brackets rather than an overbar. For L below L c rs 0.73 there 
is only a RS solution which is very close to the lattice prediction. For L greater than L c 
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FIG. 2. Plot of (R? F (L)) vs. L. The parameters are M = 1/2, g = 25, p = 1, ^ = and 
H = 0.01. The dotted line is generated by averaging over 10000 samples on a lattice of size S=40 
with A = 0.2. The error bars are found by computing the standard deviation of 10 sets of 1000 
samples. The dashed line is the RS solution, and the solid line is the RSB solution. 
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the RS and RSB solutions are different and it is clear that the RSB solution is closer to 
the lattice result. We can see that the end-to-end distance saturates at a constant value as 
L increases. This behavior is correctly predicted by the RSB solution but not by the RS 
solution. 



We now turn our attention to the quantity (Rg(L)). In Fig. |3| we plot (Rq(L)) vs. L 
using the lattice computation and also using the RS and RSB solutions. We can clearly see 
that when the RS solution differs from the RSB solution, the RSB prediction is closer to 



the lattice prediction. Again, the quantity (Rq(L)) becomes constant for large L and this 
is correctly predicted by the RSB solution. For this quantity similar results were obtained 



previously in Refs. M and [15 



It is clear from Fig. ^| and [| that the variational method with the quadratic ansatz in 
Eq. ( |4.5| ) is quite effective in describing the physical properties of the polymer chain. The 
features predicted by the lattice computation are consistent with a RS solution for chain 
lengths shorter than L c and a RSB solution for chains longer than L c . In a later section we 
will explore the physical interpretation of the variational solution and show that it is indeed 
consistent with the physics of the problem. 



B. Localized eigenstates and glassy behavior 

In this section we explore, using the lattice computation, the connection between the 
eigenstates of the Schrodinger equation and the physical properties of the polymer chain. 
We focus on the probability distribution defined as P(R, L) = Z(H, H,L)/ f Z(H, R, L)dH 
which can be interpreted as the probability of finding a closed polymer chain of length L 
which passes through the point R (for a given realization of the random potential). We 
consider this probability distribution since it gives the most direct connection between the 
chain properties and the eigenfunctions of the Schrodinger equation. In Fig. [| we plot 
P(R, L) vs. R for four different chain lengths. We also include a plot of the random 
potential sample that is used. We can see clearly that the probability distribution evolves 
from a flat distribution to distributions that are sharply peaked at various locations in the 
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FIG. 3. Plot of (Rq(L)) vs. L. The parameters are the same as those of Fig. || except that 
here we use [i = 0.3. The dotted line is generated by averaging over 10000 samples. The error bars 
are found by computing the standard deviation of 10 sets of 1000 samples. The dashed line is the 
RS solution, and the solid line is the RSB solution. 



10 20 30 40 

R 

FIG. 4. Plot of P(R,L) vs R for four values of L. The bottom most curve is the random 
potential sample that is used. From bottom to top we use L = .3, 1, 10, 20. The parameters are 
the same as those used for Fig. 
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FIG. 5. Plot of <S> 2 m (R) vs. R. From bottom up we plot the eigenfunctions with m = 0, 10, 19, 39. 
We use a lattice of size S=40 with A = 0.133. All other parameters are the same as that of Fig. ^. 

sample. Also, an important feature is that the number of peaks decreases with increasing L 
until finally there is only one peak. This implies that longer chains have a tendency to be 
found in a few favorable regions in the sample, while short chains can be found with equal 
probability almost anywhere in the sample. Another observation is that the peaks in the 
distribution, for the cases L — 1, 10, 20, are concentrated around the valleys of the random 
landscape. This shows that the chain is more likely to be found in regions of low average 
potential. As can be seen the width of the well is also important. 

The results of the previous section can be better understood by studying the properties 
of the eigenfunctions of the Schrodinger equation in a random potential. Using the eigen- 
function expansion in Eq. ( |3.7|) we can see that P(R, L) is just a sum of the eigenfunctions 
squared weighted with the Boltzman factor e -P LE ™_ g Q ^he shape of P(R, L) is dependent 
on which eigenstates have the dominant weight at a given chain length L. In Fig. [5] we 
plot a representative sample of the eigenfunctions for a typical random potential sample. In 
order to get a more general picture of the structure of the eigenfunctions we measure for 
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FIG. 6. Plot of w m vs. m. We use a lattice of size S=40 with A = 0.133 and plot w m for the 
first 75 eigenstates. We use 10000 samples and the parameters are the same as that of Fig. |2|. 

each eigenfunction the width defined by 

w 2 m = J R 2 e m {R)dR- (| R& m {R)dRj . (6.1) 

This width is closely related to the localization length £ m defined earlier. We then average 
this quantity over many realizations of the random potential and in Fig. |6] we plot w m vs. 
m. It is clear that as m is increased the width of the eigenf unctions also increase. Clearly, 
on average the ground state is the state with the smallest width. This explains why the 
the distribution P(R, L) evolves with L into fewer and fewer sharply localized peaks. This 
is just because the low energy eigenfunctions are sharply localized and they dominate the 
partition sum as L is increased. The higher excited states have typically large localization 
lengths, and these are the states that dominate the partition sum when L is short. For the 
lowest states we have verified the relation w 2 m ~ l/\E m \ . In Fig. [7| we plot w 2 m \E m \ vs. 
m for the lowest ten eigenfunctions and show that it is nearly constant in agreement with 
Eq. (pi). 



We now consider the role of the eigenfunctions on averaged quantities such as (R 2 F (L)) 



and (Rq(L)). In order to study the effect of the eigenfunctions on (R|,(L)) we truncate 
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FIG. 7. Plot oiw^lEml vs. m for the lowest ten eigenfunctions. We average over 10000 samples 
and the parameters are the same as that of Fig. ||. 

the sum over eigenstates J2o to J2o m ^ ne eigenfunction expansion in Eq. ( |3.9|) . So by 
changing the index j we can see how the average end-to-end distance depends on how many 
eigenfunctions are kept in the expansion. In Fig. |8] we plot the average end-to-end distance 
(R|>(L)) for a number of different j values. As expected the flat portion of the graph when 
L is large is entirely due to the ground state contribution. So the constant (R 2 F (L)) is just 
a consequence of ground state dominance. This is consistent with the results in Sec. Jll. 
where we showed that for long chains (R F (L)) = 2d(d+ l)(ig S ) which is indeed independent 
of L. We can also see that the curved portion of the curve corresponds to the case before 
ground state dominance when a number of tail states contribute to the partition sum. For 
small L we see that more than ten eigenfunctions are needed in order to capture the correct 
behavior of (R 2 F (L)). 

As discussed in the previous paragraphs it is clear that the evolution of the partition 
sum is dependent on the nature of the eigenfunctions. Consequently, the sample-to-sample 
variation of quantities such as (R F (L)) V and (Rq(L)) v are also crucially dependent on 
the eigenfunctions. These sample-to-sample variations are important since they provide a 
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FIG. 8. Plot of (R|(L)) vs. L. The thick solid line is computed using all the eigenfunctions. 
The thin solid line is found using only the ground state. The long dashed line is using the first five 
eigenfunctions and the short dashed line is using the first ten eigenfunctions. 

measure of the glassiness of the system. In Fig. |9] we plot (R? F (L)) V vs. L for 5 randomly 
picked samples. Notice that the spread of (R 2 F (L)) V around the average increases with 
L. In order to quantify this spread more accurately we plot the relative sample-to-sample 



fluctuations A F /(R F (L)) vs. L in Fig. [T(| It is clear again that as the the length of 
the chain is increased the sample-to-sample fluctuations also increase. Comparing Fig. [L0 
with Fig. § it is clear that the sample to sample fluctuations increase rapidly at the same 
approximate value of L at which the low energy eigenstates begin to dominate the behavior 



of (R, F (L)). The reason for this is that distributions dominated by localized states are 
strongly sample dependent. On the other hand when L is small and the high excited states 
dominate, the sample-to-sample fluctuations are small since the shape of the extended states 
do not depend strongly on the random sample. So we can conclude that for a polymer in 
short range correlated random media, the glassy characteristics are directly related to the 
dominance of localized eigenfunctions. 
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FIG. 9. Plot of (R^(L))y vs. L. The dashed lines correspond to five randomly chosen samples. 



The thick line is the average (Rp(L)) over 10000 samples. 




FIG. 10. Plot of Ap/ (Ri(L)) vs. L. The parameters are the same as in Fig. |2| and we average 
over 1000 samples. 



26 



C. Analysis of (R%(L)) 



We will now consider the end-to-end distance (R^(L)). Since one end is tethered to the 



origin it turns out that this quantity behaves very differently from the quantities (R|,(L)) 



and (Rq(L)). The replica variational approach described in Sec. |V| is not adequate in this 
case in this case and hence we will only present the results of the lattice computation. 
For short range correlations we found that the numerical method described in Sec. 



was unreliable. The reason for this is that the sum over energy eigenfunctions in Eq. ( |3.3|) 
is unstable since the overlap <& m (R)<$> m (R') (for short range correlations) is a number on 
the order of exp(— \R — R'\/£ m ) which is typically extremely small. However, we were able 



evaluate Eq. (<T3) accurately for small £ by solving the Schrodinger equation on a lattice 
using a fourth order Runge-Kutta algorithm with a very small time step (t ~ 10 -3 ). We 
performed the computation on a large system (S = 160, iV = 800) in order to minimize 
finite size effects. 



In Fig. |TT]we plot (R|,(L)) vs. L and also include (Rr(I))y vs. L for five typical samples. 
Notice that for a given sample the end-to-end distance (Ry(L))y has a tendency to change 
rapidly and then remain constant. As discussed in Sec. |I| this is simply a consequence 
of hopping between localized states. That is, the free end of the polymer rapidly finds a 
deep well (which is also a localization center) and stays there. Upon averaging over 10000 



random samples we find that (Ry(L)} grows linearly with L for small L but scales faster 
than diffusion for large L. The diffusive behavior when L is small is not surprising since a 
short chain should not be affected by the random media. The claim based on Flory type 



arguments [§]|| is that (R^(L)) ~ L 2 / (log(L)) 7 when L is large. However, we were not able 
to extract a consistent exponent 7 for the range of L that we explored. The reason for this is 
that our data is not accurate enough to detect the precise logarithmic correction to scaling. 



However we went ahead and tried to find the best power law fit (R^(L)) oc L u for large L. 
We found that for 5 < L < 25 the exponent v « 1.74 yielded an excellent fit to the data. 
This result is consistent with the sub-ballistic prediction as it is very close but slightly less 
than ballistic scaling (L 2 ). 
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FIG. 11. Plot of (Rji(L)) vs. L. The full circles are generated by averaging over 10000 samples, 
and the error bars are found by computing the standard deviation of 10 sets of 1000 samples. The 
dashed lines are plots of (Ry(L))y vs. L for five typical samples. The parameters are the same as 
in Fig. ||, only that here the system size is four times larger (S = 160) and we set n = 0. We use a 
lattice of 800 sites with A = 0.2. 
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FIG. 12. Plot of x c vs. L. The parameters are the same as those used in Fig. 0. 



D. Analysis of the variational solution 



So far we have seen that the glassy characteristics of a polymer in a short range correlated 
random potential is closely related to the dominance of low energy eigenfunctions. We also 
know that the variational solution possesses a RS solution for L < L c and an RSB solution 
for L > L c . It is well known that replica symmetry breaking is typically associated with 
glassy behavior, and for our model we show that the onset of RSB is precisely when the 
system begins to exhibit glassy behavior. The variational parameter that best reveals the 
transition between RS and RSB is the break point x c . If x c = 1 then there is only an RS 
solution, and if < x c < 1 then that corresponds to an RSB solution. In Fig. [12] we plot 
x c vs. L using the same parameters that were used in Fig. 0. We can see that onset of the 



RSB solution is at L c pa .73. If we compare this result to the plot of Ap/ (R|,(L)) vs L in 



Fig. [lCj, we see that near L c pa 0.5 the sample-to-sample fluctuations begin to rise rapidly. 
This result provides strong evidence that when the RSB solution is valid the polymer chain 
does indeed exhibit glassy behavior. 

In Ref. approximate analytic solutions to the variational equations were found. This 
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FIG. 13. Plot of A vs. L. The parameters are the same as those used in Fig. |2[ 

was for the case of large L, small /x, and a delta correlated random potential (£ — > 0). It 
was found that 



L V(2vr) 



J(Z-2 \ -V(4-d) 

" ^ +4 M d |log/i| d - 2 



(6.2) 



We checked the L dependence of x c in Fig. [12| and indeed we found that x c oc 1/L. In fact, 
we found that this 1/L dependence is quite robust as it holds whenever there is a 1-st ep 
RSB solution. In a later section we will analyze the physical consequences of this behavior 
and show that it can be simply explained via the path integral mapping to the Schrodinger 
equation. 

In Fig. |T3] we plot the parameter A vs. L. The discontinuity is of course at L c , after 
which we plot only the 1-step RSB solution. Notice that for L > L c A is essentially constant, 
which by Eq. 



implies that (R, 2 F (L)) is also constant. This result is consistent with the 
approximate analytic solution in Ref. 0, where it was found that 



A 



d 4/(4-d) 



(27r)2d/(4-d) 

which is a constant independent of L. 



-^M)^y^ d \g\\og^- d \ 



(6.3) 
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VII. PHYSICAL INTERPRETATION OF THE 1-STEP RSB SOLUTION 



In this section we study the physical interpretation of the replica symmetry breaking 
solution. Our purpose is to see if the underlying physical picture predicted by 1-step RSB 
is indeed consistent with the presence of exponentially localized eigenstates. The following 
analysis is valid for a very long polymer (large L) when the system becomes glassy. We 
begin by evaluating the replicated partition sum defined as 

"Ra(L)=Ro n 

'R a (0)=R a a=1 

where h n is the quadratic trial Hamiltonian in Eq. (4J5). Since h n is quadratic the path 



pRa(L) = R a n 

Z n ({R a }) = Y[[dR a }exp(-(3h n ), (7.1) 



integrals can be evaluated analytically and the final result can we written in the form 

Z n ({R a }) = const, x exp ^ Q^R a ■ R^j . (7.2) 

The details of this calculation along with the relationship between the matrices Q a t, and p a b 
are given in the Appendix. Now, since p a b was parametrized according to the 1-step RSB 
scheme, it implies that Q a b can also be parametrized in the same way. 

Mezard and Parisi |20| discuss the interpretation of a representation of the form ([7.2|) 
for the case of directed polymers. In particular they show how to deduce the structure of 
the probability distribution 



P V (R) = Z V (R,L)/ J dRZ v (R,L), 



(7.3) 



which is the probability of finding a polymer loop that passes through R for a given re- 
alization (which we denote by V) of the the random potential. Here Zy(R,L) is just the 
partition sum Z(R, R; L) as given in Eq. (|2.4j ). This probability is related to the replicated 
partition function given in Eq. f[7.1|) by 



P V {R) = lim / dR 2 ■ ■ ■ dR n Z n ({R a }) Kl=R . (7.4) 

n^O J 

Mezard and Parisi's analysis has to be adapted for the case of real (non-directed) polymers 
of length L in a random potential which is independent of time. The changes will be pointed 
out below. 
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If Q a b is parametrized by {q, q(x)} such that 

{go x < x c 
(7.5) 
qi x > x c , 

one proceeds to obtain Py(R) by the following procedure: 

1. For each sample ( a realization of the random potential) generate a random variable R 
which is picked from the distribution 

2. Consider a set of "states" labeled by the index a whose physical meaning will be 
elucidated shortly. Each of these states is characterized by a weight W a and a position 
variable R a . Given R , the variables R a are an infinite set of uncorrelated random 
variables distributed according to 

V(Ru R 2 , ) = II To 7 1 ^ ex P f- ( ^~ R °D • (7-7) 
The distribution of weights will be discussed below. 

3. Given these "states" for a given sample, The probability distribution Py(R) for that 
sample has the form 

The weights W a are given in terms of some "free energy" variables f a : 

Wa = ^tfM (7 . 9) 
L 7 exp(-/3/ 7 ) 

These free energy variables are chosen from an exponential distribution 

P\f a ]<xexp(xcPf a )0(f-f), (7.10) 

where / is an upper cutoff. 
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What is the meaning of these variables in the present case? To determine the weights 
W a we compare Eq. ( |7.8| ) to the eigenfunction expansion given in Eq. (|3. 7[). From Eq. ( p.7| ) 
together with Eq. ( [7.3| ) it becomes clear that 

P y (R) = ^A Q |$ Q (R)| 2 , (7.11) 

a 

where 

exp(-pLE a ) 
E 7 exp(-/3LE 7 ) 

Comparing Eq. ( |7.8| ) and Eq. ( |7.12| ) it becomes obvious that 

W a = A a , (7.13) 

and 

^(R)ocexp f- ^r^r V ( 7 - 14 ) 



Hence, the "states" labeled by a are in our case the actual eigenstates of the imaginary time 
Schrodinger equation ( |3.4| ). These are localized tail states centered at position H a with an 
associated "weight" W a . Thus the 1-step RSB solution approximates the tail states by a 
fixed Gaussian form. The width of these Gaussians (w ) as defined by Eq. ( |6.1| ) satisfies 
Wq = d(q — qi), which for large L can be shown to converge to Wq ~ d/(2pT^/\M) (see 
Appendix). Since A becomes constant for large L then so does the width u> . Thus the 
1-step RSB solution approximates all the localized states by a Gaussian of constant width 
Wq, which is the typical size of a Gaussian chain embedded in a random potential. It also 
becomes evident that the free energies f a are equal to LE a . This make sense if we think 
of \E a \ as representing the binding energy per monomer, and thus f a = LE a represent the 
total energy of the chain. Alternatively one can start with the quantum particle picture and 
identify f a with E a but then in the transformation from a quantum particle to a polymer 
one has to replace (3 by PL, and thus one ends up with the same expressions. 

These arguments lead us to expect that within the 1-step RSB scheme, the energy vari- 
ables E a are independent random variables taken from an exponential distribution: 

P[E a ] oc e pLxcEa 6(E - E a ) , (7.15) 
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with E being some energy scale determined by the upper cutoff of the tail region. We will 
now argue that the distribution given above is just the expected distribution of ground-state 
energies i.e. the probability of finding the lowest energy level to have energy E. We first 
review some very basic results of extreme value statistics as presented in Ref. p4 |. Given 



K independent and identically distributed random variables E iy pulled from a distribution 
of the form 

P(E) = j^ex V (-B\E\ s ), (7.16) 

the probability that the lowest of the K energies is E (for E — > — oo and K — > oo) is given 
by 

P(E) oc exp [B5\E C \ 5 - 1 E] , (7.17) 

where 

^-fffly" (7,8) 



B 

The value of E c , the lowest energy expected to be attained in K trials, is easily obtained 
from 

dEP(E) ~ l/K. (7.19) 

The reason why we chose a distribution of the form given in Eq. (|7.16|) is that in d = 1 the 
probability P(E) is known to have that form exactly for the case of delta correlated random 
potentials (see fTSfl). For d > 1 Lifshits ^5] argued that the form given by Eq. ( |7.16| ) is also 
valid. Our goal now is to see if the distribution Eq. ( |7.15| ) derived using the 1-step RSB 
solution is indeed consistent with the distribution Eq. ( |7.17| ) predicted using extreme value 
statistics. 

Comparing Eq. ( |7.15| ) and Eq. ( [7.17| ) we find that for consistency the break point should 
satisfy 

x c = -^B^ilog^K))^ 5 . (7.20) 
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Notice that the 1/L behavior of x c is exactly the same as was found analytically for large L 
in Ref. [7J] and numerically for any L > L c in the present work. 

We can go further by using the fact that the number of energy levels K, within a 
fixed energy interval is directly proportional to the system size, which in our formulation is 
effectively determined by fi. Assuming log(K) oc | log(/i)| and comparing to the approximate 
solution for x c in Eq. (|6.2|) we find that 5 = (4 — d)/2 and B oc 1/g. Now P{E) is just 
proportional to the density of states p(E), and it is known exactly in one dimension. Indeed, 
when d = 1, S = 3/2 and B oc 1/g. For 2 < d < 4, 5 agrees with the result derived by 



Lifshits ||25|| . Hence, the exponent 6 and the disorder dependence of B is correctly predicted 
by the 1-step RSB solution. The exponent a cannot be found since the statistics of the 
lowest energy depends only on the exponential tail. 

The 1-step RSB solution is thus in agreement with Eq. ( |7.17|) . Since x c oc 1/L then the 
distribution of the energies predicted from Eq. ( |7.15|) is independent of L. Note also that the 
value of E c which gives a typical lowest energy for a given volume of the system coincides 
with the expected value given in Eq. ( |2.7] ) when one substitutes for B ~ 1/g, S — (4 — d)/2 
and In K ~ In V in Eq. ( |7.18| ). The localization length which is the width of a localized state 
is known to be related to the energy of the state by 

e ~W^W- (7 ' 21) 

If we use E c for the value of the lowest energy in a volume V we get an estimate for Rp, the 
radius of gyration of a free chain which is in agreement with Eq. (|2.6| ) . 

When L — > oo most of the weights W a — > and only the weight with the lowest E a 
will contribute significantly to the distribution Py(FL). So in this limit Py(R) will be 
approximated by a single Gaussian located somewhere in the sample. This result is consistent 
with ground-state dominance as predicted by the corresponding eigenfunction expansion. 

We now consider the distribution P(R a ) given in Eq. ( ff77|) . This is just the distribution 
for the localization centers R Q for a given value of Ro. Hence, we can calculate the average 
distance between the localized states for a given sample. We find that the width w of the 
Gaussian P(R a ) satisfies w 2 = d(qi — go)- For small /i and large L we find (see Appendix) 
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that w 2 ~ d/((3fiLx c ). Using the analytic approximation for x c we find that 

1 



W - ~ _^2/(4- rf )| l og ( /x )|(^2)/(4- C? )_ (? _ 22) 



,2 



We can also find the average distance from the origin to a localization center R Q . To do 
this we have to compute the average probability of find a localization center at R a . This is 

P(R a ) ex J rfR exp (-R*/(2g )) exp (-(R a - R )7(2(<?i - go)) 

= exp (-R^/(2 gi )) . (7.23) 

If we estimate the average distance to a localization center to be the width of the Gaussian 
P(R a ), then w 2 = dgi « d/(j3[j,Lx c ) which is the same as the average distance w between 
tail states derived above. This is because the fluctuations of Ro are small compared to the 
fluctuations of H Q — Rq when the volume is large. Notice that this result is consistent with 



our calculation for (Rq(L)) in Eq. ( |4.1U| ), since for L large (Rq(L)) ~ d/(j3/iLx c ). These 



quantities should be close since both give an approximation to the average distance to the 
ground state localization center. 

So far we have studied the physical interpretation of the 1-step RSB solution, but recall 
that for L < L c there was only a RS solution. It is interesting to study the physical 
implications of this RS solution. In this case the probability distribution -Py(R) will have 
the simple form 

P ^ R ) = ,<> r~ 1 sm exp f- ( ^~ R °0 > ( 7 - 24 ) 

(2?r(g - g )r /2 V 2 (<? ~ 9b) / 
where the variable Ro is taken from the distribution 

p < r »^ (-!) ■ < 7 - 25 » 

Notice that for a given sample Py(R) is approximated by a Gaussian shifted from the origin 
by a distance Ro- The width of this Gaussian can be shown to be essentially independent of 
disorder in the limit /i — > 0, as the width wq satisfies Wq ~ d^/SfiL)' 1 , which is the result for 
the case of zero disorder. The average shift from the origin can be estimated by the width 
wi of the Gaussian in Eq. ( |7.25| ), which satisfies wf ~ g(3^ d+2 ^ 2 L^ d+2 ^ 2 ^■ d ~ T) l 2 , for the case 
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of delta correlations and large volume (// — > 0). This result is similar to that found for the 
case of long range correlated random media in which it was found that the only effect of 
disorder is to shift the center of mass of the polymer chain from the origin. 

We can conclude that the physical interpretation of the 1-step RSB solution is consistent 
with that given by the eigenfunction expansion. Clearly, the 1-step RSB solution captures 
the localized states and also correctly predicts some important features of the eigenvalue 
distribution. However, there are differences and these reveal the limitations of the 1-step 
RSB solution. For example, all the localized states are approximated by the same Gaussian 
profile when in fact the localization lengths should increase with energy. 

VIII. CONCLUDING REMARKS 

The findings in this paper allow us to give a fairly complete description of a polymer 
in random medium. We find that the size of a short polymer chain behaves as if there is 
no disorder, but that the position of the center of mass of the chain is strongly affected 
by the random media. Applying the replica variational method we find that short chains 
are well described by a RS solution. Also, using an eigenfunction expansion we show that 
the partition sum of short chains is dominated by extended eigenstates. For long chains 
the physical picture becomes more interesting and complex. We find that long chains are 
likely to be localized in regions of the sample where there is very low potential energy. These 
regions correspond to localized tail states of the corresponding Schrodinger equation. As the 
length of the chain increases the number of dominant conformations decreases until finally 
the chain is essentially localized in one small region, consistent with ground state dominance. 
In terms of the mapping from a chain to a quantum particle, a big chain length corresponds 
to a very low temperature for the quantum particle. 

We show that the onset of the variational RSB solution to the stationarity equations 
corresponds to the dominance of localized tail states. This 1-step RSB solution correctly 
describes the glassy characteristics of the polymer chain such as large sample-to-sample 
variations. We have demonstrated the clear physical picture associated with the abstract 
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RSB solution. 

Our analysis also suggests a direct connection between the size of a polymer in a random 
potential with short range correlations, and the localization length associated with the lowest 
energy state in a system of finite volume V. Given the density of states p(E) we can estimate 
the lowest energy E c by 



dEp(E) ~ 1/V. (8.1) 

The size of a free Gaussian chain, in the long chain limit, is then given by 

R F ~ \E C \- 1 ' 2 . (8.2) 

It is remarkable that starting from an annealed average (the density of states involves the 
average of the partition function) and using extreme value statistics plus the known relation- 
ship between the localization length and the value of the tail state energy, we can predict 
the behavior of a quenched average like the chain size. 

It is well known that glassy systems are notoriously difficult to simulate using the Monte 
Carlo method since equilibration times are exceedingly large. In the context of a polymer 
chain in a random medium, we expect that a long chain will get trapped in the deep wells 
of the random potential and any local updating procedure, such as the Metropolis scheme, 
will not be able to find the true minimum energy configurations in a reasonable time. A 
workaround is to use a simulated annealing procedure Of course, these traps are just 



those regions of the sample where an eigenstate would be localized. They correspond to 
localized states with energies above the true ground state that exists for a finite-size system. 
Since the partition sum for a long chain is dominated by localized tail states those chain 
configurations that are in the vicinity of such a state are overwhelmingly more favorable than 
other configurations. All the classic characteristics of a glassy system emerge (like trapping, 
metastability, aging), and they have a clear physical origin in terms of the localization 
picture. 
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APPENDIX: 



In this appendix we derive the relation between the matrix p ab (see Eq. (|4.7| ) ) char- 
acterizing the variational hamiltonian and the matrix Q a f, defined in Eq. ( |7.2p . We start 
from 



Z n ({R' a },{R a };L) 



R a (L)=R a n 

exp(— (3h n ), 

R a (0)=R^ a=1 



(Al) 



with 



kn = lJ M ( 



'd R a (u) 
du 



+ XRliu) 



+ 



— / du f du'y2pabR a (u) ■ R b (u'). 



(A2) 



R a (L)=R Q ri 

, Y[[dR a (u)}exp(-{3h n ({\ c })), (A3) 

Z I -/R a (0)=R' a a=1 



Let Ooh be an orthogonal matrix which diagonalizes p a b by a similarity transformation, and 
let p c , c — 1, ■ ■ ■ ,n be the eigenvalues of p a b. We can express Z n in the form 

Zn = / n exp 

c 

with 

hn({K}) = \£ duJ2 M (^7^) + AR*(lt) + J2 V-Pc/mOac K " Ra(u) . (A4) 

The path integral can now be carried out using Eqns. (3.39-3.41) in [IB|. The result is 



] Jr/A r cxp - 



(A5) 



where 



$ = A(^R^ + ^Rl 2 )+fi(^R a + ^Rl).V a + 2C^R a -Rl + J D^V^ (A6) 
with 

A = 1 coth (l VA/m) ( A7) 



B = VXM [cosh (Ly/X/M) - ll [sinh (l^X/M 



-i 



C = —VXM [(sinh [Ly/\/M 
D = B — LA/2, 

V a = i ^ v/-j9 c /(/3L)0 QC A c . 



(A8) 
(A9) 
(A10) 
(All) 
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Notice that 

a a 

We now put R a = R4 and integrate out over A c to obtain 

Z n {{R a }, {Rj; L) = Afexp ^-2/3 ^[{A + C)5 ab - [3B 2 M ab ]R a ■ R b j , (A13) 



with 



If X 2 L " ' 



M = [ 1 - — p _1 . (A14) 

Here we denoted by 1 the unit n x n matrix, and p stands for the matrix p ab . Comparing 
with Eq. ([FJ) we find that 

Q ab = jg{[(A + C)l - PB'M]- 1 }^. (A15) 

In the Appendix of Ref. [[2(J one can find a formula for the inverse of a hierarchical matrix. 
For the one-step RSB case, the inverse of a matrix p = {p,Po,Pi} is given by q = p 1 where 

«=r l -n rA 1 —- ^ _ (pi - po)(i - Zc) \ (A16) 

(P -p X c -Pl{l ~ X c )) V P -P0X c -Pl{l - Xc) P-Pl J 

Qo = -jz w^, (A17) 

[p-pQXc-pi.il - X C )Y 

?i = - ~ t: r + ~ t: r • ( Al8 ) 

p-_p ^c-Pi(l-^c) \P-Pi p-p x c -pi{l-x c )J 



In the limit of large L we have 



A^^v/AM (A19) 



B -> VAM (A20) 
C -> (A21) 



LA 



D -> + VAM. (A22) 

2 

In addition we use the one step breaking results found in Ref. [0] . Evaluating the matrix Q 
in the limit of small /x we find to leading order 
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? = T^TT + 77TT + l -F= + ■■■ ( A23 ) 

Lx c Pfi L(3/j, 2 2[3Vm\ 

* = I | ? + - (A24) 



" = + W + (A25) 



with the parameters A, and x c given in Eq. (|4.13[) and ( |4.14|) respectively and so is given 
approximately by [^] 

so = const, x g<?-WpL(j3VM)^W f i d ' 2+1 |ln/i| {d+2)/4 . (A26) 

We have also used the fact that £ = A in this limit. 
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